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ABSTRACT 

Transits of hot Jupiters in X-rays and the ultraviolet have been shown to be both deeper and more 
variable than the corresponding optical transits. This variability has been attributed to hot Jupiters 
having extended atmospheres at these wavelengths. Using resolved images of the Sun from NASA’s 
Solar Dynamics Observatory spanning 3.5 years of Solar Cycle 24 we simulate transit light curves of a 
hot Jupiter to investigate the impact of Solar like activity on our ability to reliably recover properties 
of the planet’s atmosphere in soft X-rays (94 A), the UV (131-1700 A), and the optical (4500 A). We 
find that for stars with similar activity levels to the Sun, the impact of stellar activity results in the 
derived radius of the planet in soft X-ray/EUV to be underestimated by up-to 25% or overestimated 
by up-to 50% depending on whether the planet occults active regions. We also find that in up-to 
70% of the X-ray light curves the planet transits over bright star spots. In the far ultraviolet (1600 
& 1700 A), we find the mean recovered value of Rp/R^ to be over-estimated by up-to 20%. Eor 
optical transits we are able to consistently recover the correct planetary radius. We also address the 
implications of our results for transits of WASP-12b and HD 189733b at short wavelengths. 

Subject headings: stars: activity - stars: coronae - planets and satellites: individual HD 189733b, 
WASP-12b) - planets and satellites: atmospheres - stars: starspots 


1. INTRODUCTION 

Understanding the impact of stellar activity on transit 
and radial velocity data is not only fundamental for find¬ 
ing Earth-sized planets but also in determining the com¬ 
position of both the largest and smallest exoplanets. Eor 
instance, by observing transits at multiple wavelengths 
the composition of the atmosphere of the planet can be 
inferred through transmission spectroscopy (see for ex¬ 
ample Sing et al. 2011; Nikolov et al. 2013; Stevenson 
et al. 2014; Gibson et al. 2012). This process requires 
an accurate and precise measurement of the planetary 
to stellar radius ratio (Rp/R*) at various wavelengths. 
It is therefore crucial that the effects of stellar activity 
on transit light curves are well understood since activity 
will alter the shape of the transit light curve. 

The impact of stellar activity on transit light curves 
at different wavelengths has been used to not only learn 
about the planet but also the host star itself. At opti¬ 
cal wavelengths, bumps in the transit light curve caused 
by the planet occulting star spots can be used to infer 
the presence of magnetically active regions on the stel¬ 
lar surface. By tracking the locations of these bumps 
in consecutive transits the spin-orbit-alignment between 
the orbital plane of the planet and rotation axis of the 
star can be derived and this value is in good agreement 
with Rossiter-Mclaughlin measurements (Sanchis-Ojeda 
& Winn 2011; Winn et al. 2010). Eor misaligned sys¬ 
tems it may be possible to track changes in the location 
of bumps in the light curve over time to recover stel¬ 
lar butterfly patterns to investigate whether other stars 
exhibit such cyclic evolution in the distribution of star 
spots (Llama et al. 2012; Sanchis-Ojeda et al. 2013). Un¬ 
derstanding the properties of stellar butterfly patterns 
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will provide insight into the stellar dynamo (Berdyugina 
2005). 

Asymmetries in the near-UV (2800A) light curve of 
WASP- 12b observed using Hubble Space Telescope (HST) 
have revealed an early-ingress and also a deeper transit 
when compared to the optical light curve (Eossati et al. 
2010; Haswell et al. 2012). However, the end of the near- 
UV transit coincided with the optical egress implying the 
presence of additional, asymmetrically distributed ab¬ 
sorbing material in the exosphere of the planet in the 
near-UV. 

One potential explanation for a deeper, asymmetri¬ 
cal transit is that heavily inflated hot Jupiters such as 
WASP-12b can overflow their Roche lobes (Gu et al. 
2003; Li et al. 2010; Ibgui et al. 2010) which may result 
in an optically thick accretion stream from the planet 
onto the star (Lai et al. 2010). Another solution is that 
the early ingress is the signature of the stellar wind im¬ 
pacting on the magnetosphere of the planet creating a 
magnetospheric bow shock (Vidotto et al. 2010). This 
idea was modeled using Monte-Carlo radiative simula¬ 
tions in Llama et al. (2011) where they were able to fit 
the near-UV observations of Eossati et al. (2010) with 
a bow shock model that suggests the planet’s magneto¬ 
sphere extends out to a distance of ^ 5Rp. 

Hot Jupiters orbiting much closer to their parent star 
are exposed to higher levels of stellar irradiance than the 
gas giants in our Solar System. This increase in stellar 
irradiance may heat the exosphere of the planet lead¬ 
ing to mass loss (see for example Lammer et al. 2003; 
Lecavelier Des Etangs et al. 2004). It may be possible to 
detect such signatures in temporal variations of transit 
light curve depth and duration. 

Observations of HD 189733b at multiple wavelengths 
have been taken searching for evidence of an inflated 
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TABLE 1 

Wavelengths observed by the Atmospheric Imaging Assembly (AIA) on-board NASA’s Solar Dynamics Observatory (SDO). 

Adapted from http://aia.lmsal.com/public/instrument.htm. 


Wavelength 

(A) 

Regime 

Primary ions 

Source region 

Temperature (K) 

Limb 

bright / dark 

94 

Soft X-ray / EUV 

Fe XVHI 

Flaring and quiet corona 

6.3 X 10^ 

Brightened 

131 

EUV 

Fe VHI, XX, XXH 

Flaring and quiet corona 

4 X 10^, 10'^, 1.6 X 10'^ 

Brightened 

171 

EUV 

Fe IX 

Quiet corona and upper transition region 

6.3 X 10^ 

Brightened 

193 

EUV 

Fe XH, XXIV 

Corona and hot flare plasma 

1.2 X 10®, 2 X 10'^ 

Brightened 

211 

EUV 

Fe XIV 

Active regions of the corona 

2 X 10® 

Brightened 

304 

EUV 

He H 

Chromosphere and transition region 

5 X 10^ 

Brightened 

335 

EUV 

Fe XVI 

Active regions of the corona 

2.5 X 10® 

Brightened 

1600 

EUV 

C IV, continuum 

Transition region and upper photosphere 

10®, 5x10® 

Darkened 

1700 

EUV 

continuum 

Temperature minimum, photosphere 

5 X 10® 

Darkened 

4500 

Optical 

continuum 

Photosphere 

5 X 10® 

Darkened 


planetary atmosphere and mass loss from the planet. 
Ben-Jaffel & Ballester (2013) reported a potential early 
ingress and also a deeper transit depth of the C ll emis¬ 
sion line at 1335 A which they attribute to the mag- 
netospheric bow shock occulting star light. Their MHD 
modeling predicted the bow shock will have a stand-off 
distance of ^ 17Rp. 

Lyman-(a transits have been used to study both atmo¬ 
spheric escape and inflated planetary atmospheres on a 
number of hot-Jupiters. Vidal-Madjar et al. (2003) ob¬ 
served transits in Lyman-a (1215 A) of HD 209458b and 
found the first evidence of an extended exoplanetary at¬ 
mosphere with a Lyman-(a transit depth ^ lOx deeper 
than the corresponding optical depth. Linsky et al. 
(2010) also observed HD 209458b and found a deeper 
depth in their C II and Si ill transit light curves, sug¬ 
gesting mass loss from the planet. Lecavelier Des Etangs 
et al. (2012) observed transits of HD 189733b using 
HST/STIS in Lyman-a at two epochs. In 2011 Septem¬ 
ber they found an additional 14.4% absorption and also 
hints of an early ingress when compared to their obser¬ 
vations in 2010 April, which were similar to the optical 
light curve. They attribute this variability to the stellar 
wind directly influencing the evaporation rate of the up¬ 
per atmosphere of the planet. Observations of 55 Cancri 
b have shown a partial eclipse in Lyman-a that coincides 
with the predicted transit time, suggesting the planet 
has an extended atmosphere that grazes the stellar disk 
(Ehrenreich et al. 2012). Kulow et al. (2014) observed the 
Hot-Neptune GJ 436b in Lyman-a and found not only 
a deeper transit, but also that the deepest part of the 
Lyman-(a transit occurred after the egress of the optical 
transit, suggesting this planet has an extended Hydrogen 
tail. 

Llama et al. (2013) investigated the impact of the vary¬ 
ing stellar wind conditions experienced by HD 189733b 
as it orbits around its parent star by modeling the prop¬ 
erties of the stellar wind from the surface of the star to 
the orbital distance of the planet using observed mag¬ 
netic maps of the star obtained by Eares et al. (2010). 
They found that variability in the density and pressure 
of the stellar wind will alter the orientation and density 
of the resultant bow shock, which in turn caused the pre¬ 
dicted near-UV light curves of HD 189733b to be highly 
variable in both depth and duration. 

Poppenhaeger et al. (2013) observed multiple X-ray 
(lA - 63A) transits of HD 189733b with Chandra and 
found a combined transit depth of 6 - 8%, which is ^ 
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Fig. I. — Top: Total Solar Irradiance (TSI) data from NASA’s 
VIRGO instrument on the SOHO instrument for the past ten 
years, covering Solar cycle 23 and 24 (Frohlich & Lean 1988). Bot¬ 
tom: The Solar butterfly pattern showing the fraction of the Solar 
disk covered in sunspots, produced using the data from Hathaway 
(2010). The blue region shows the timespan for this investigation. 
During this time the Sun was moving from activity minimum to¬ 
ward maximum, and the active belts have migrated from ±35° to 
±5°. 

3x deeper than the depth in the optical of 2.4%. They 
attribute this additional absorption to the planet having 
an extended atmosphere that is opaque to X-rays but 
transparent to optical light. They note however, that 
given the difflculty of observing transits in X-rays it is 
not possible to rule out the increased transit depth being 
caused by stellar activity without further observations. 
Eigure 9 of Poppenhaeger et al. (2013) shows how the 
derived transit depth from all the possible combinations 
of their X-ray light curves varies between 2.3% - 9.4%. 

In this paper we use resolved images of the Sun in 
soft X-rays, the extreme ultraviolet (EUV), the far ul¬ 
traviolet (EUV), and the optical to simulate transit light 
curves to investigate the extent to which stellar activity 
impacts light curves at these wavelengths and to deter¬ 
mine how reliably we can recover properties of a transit¬ 
ing hot Jupiter. 


2. SOLAR OBSERVATIONS AND TRANSIT MODEL 

Observations of the Sun present a unique opportunity 
to investigate the effects of stellar activity on exoplanet 
transit light curves. NASA’s Solar Dynamics Obser- 
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Fig. 2. — Top: Two sets of simultaneous images of the Sun obtained from NASA’s SDO spacecraft. Over plotted on each image is the 
trajectory of a simulated hot Jupiter with Rp/R* = 0.1 and b = —0.3. The left panel shows images where there was little activity along 
the transit chord, whilst the right panel is an example of when the transit chord intersects active regions. Bottom: The resultant transit 
light curves are color coded by wavelength. At soft X-ray/EUV wavelengths the Sun is limb-brightened and so the transit is much deeper 
at ingress and egress than mid-transit. The optical transits that intersect active regions (right) show a bump in the light curve. At shorter 
wavelengths where active regions appear bright, the light curves exhibit a dip as the planet occults an active region. 


vatory (SDO) has been observing the Sun since early 
2010. The Atmospheric Imaging Assembly (AIA) on¬ 
board SDO images the solar disk in ten wavelengths ev¬ 
ery ten seconds spanning the soft X-ray/EUV through 
to the optical (with a spatial resolution 1" and a CCD 
of 4096 px X 4096 px). A summary of the wavelengths 
observed by AIA with their emission is listed in Table 1. 
We note that it is not currently possible to observe stars 
other than the Sun in the EUV; however, we include the 
findings of these wavelengths for completeness. 

Using the Virtual Solar Observatory IDL routines pro¬ 
vided by SolarSoft^, we downloaded AIA data of the Sun 
at all ten wavelengths. We acquired images at a cadence 
of one image every 24 hours between 2011 January 01 
and 2014 September 01. In total we obtained ^ 1300 
images at each of the ten wavelengths. We checked each 
image for saturated pixels using the NSATPIX keyword 
in the EITS header of each image and found the frac¬ 
tion of saturated pixels in the images to be negligible. 
The most heavily affected image had a saturated area 
< 0.03% meaning that saturation will not affect our sim¬ 
ulated light curves. 

The top panel of Eigure 1 shows the composite To¬ 
tal Solar Irradiance (TSI) between 2004 - 2014 (Erdhlich 
& Lean 1988; Erohlich 2000). The bottom panel shows 
the Solar butterfly pattern during this time (Hathaway 

^ http://www.Imsal.com/solarsoft/ 


2010). The shaded blue region identifies the timespan of 
this investigation and demonstrates that between 2011 
January and 2014 September the Sun has moved from 
activity minimum toward activity maximum and the ac¬ 
tive belts migrated from ±35° toward the equator. 

The top panel of Eigure 2 shows simultaneous images 
of the Solar disk from AIA in eight wavelengths, namely, 
soft X-ray/EUV (94 A), EUV (131, 171, 193, 211, & 
304 A), the EUV (1700 A), and the optical (4500 A) on 
two separate dates (left and right panels)^. The optical 
and EUV images demonstrate how the Sun appears rela¬ 
tively inactive, with only a few active regions present on 
the surface. At short wavelengths the Sun appears more 
active than when viewed in the optical. The EUV/EUV 
images reveal how relatively small sun spots in the op¬ 
tical correspond to very bright, large, extended regions 
at short wavelengths. It is this increase in activity and 
active region coverage at shorter wavelengths that will 
have the greatest effect on exoplanet transit light curves. 

We simulated the transit light curve of an exoplanet 
over each image of the Solar disk. To compute the light 
curves we adapted the transit code developed in Llama 
et al. (2011) and Llama et al. (2013). This code was orig¬ 
inally created to simulate light curves of a planet tran¬ 
siting over a limb-darkened stellar disk. Rather than 

^ For brevity, we do not show 335 A and 1600 A as these results 
are similar to 304 A and 1700 A respectively. 
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simulating a stellar disk, we were able to use the images 
from SDO. We assumed the planet to be a hot Jupiter 
with Rp/R^ = 0.1. To enable a comparison between the 
light curves at each wavelength we assumed the radius 
of the planet to be the same in all the simulations, i.e., 
the size of the planet’s atmosphere does not change with 
wavelength. 

On the Sun, the activity belts are confined to within 
±35° of the equator. We therefore chose two values for 
the impact parameter, b. Firstly, 5 = 0, so that the 
planet transits across the equator of the star and the 
transit chord will not intersect the activity belts and 
b = —0.3 so the planet occasionally transits over active 
regions. 

Finally, we assumed the planet to be completely dark, 
i.e., any region of the stellar disk occulted by the planet 
will have a flux of zero. We assumed a high cadence tran¬ 
sit observation such that each light curve is comprised of 
one hundred data points, and the semi major axis of the 
planet to be a ^ 0.02 an so that the planet transits once 
every 24 hours. As the duration of the transit event is 
much shorter than the rotation period of the Sun, we 
assumed the Sun to be static in each of our simulations. 
At each step, the location of the planet over the stellar 
disk was computed and the loss in flux caused by the 
planet occulting star light was calculated to produce the 
simulated light curve. The effects of lower (and more 
realistic) cadence observations are discussed in Section 
4. 

3. RESULTS 

3.1. Impact of stellar activity as a function of 
wavelength 

The top panel of Figure 2 shows sample transit paths 
of a hot Jupiter with Rp/R^ = 0.1 and b = —0.3 across 
the Solar disk in multiple wavelengths. The images reveal 
that in the soft X-ray and EUV regime, the stellar disk is 
limb-brightened rather than limb-darkened. This differ¬ 
ence in the distribution of star light will change the shape 
of the transit light curve. The final column of Table 1 
lists whether the Sun is limb-brightened or limb-darkened 
for each of the observed wavelengths. The images in the 
left panel of Figure 2 are examples of when there is little 
activity along the transit path. As such, the simulated 
transit light curves will be examples of nearly quiescent 
transits. The images in the right panel of Figure 2 show 
active regions on the disk of the Sun intersecting the 
transit path causing the transit light curves to deviate 
from the nearly quiescent shape. 

The bottom panel of Figure 2 shows the correspond¬ 
ing simulated transit light curves (color coded by wave¬ 
length). The left panel shows nearly quiescent light 
curves where the planet has not transited over large ac¬ 
tive regions. In these high energy transits, a dip during 
both ingress and egress is registered in the light curve. 
The transit is also shallower at mid-transit when com¬ 
pared to the optical light curve. This is a consequence 
of the stellar disk being limb-brightened at these wave¬ 
lengths, so the planet occults a much larger proportion 
of the total star light as it transits over the limbs of the 
star, as opposed to the optical case where the light is 
more evenly distributed over the entire disk of the star. 

At short wavelengths, active regions have a higher con¬ 


trast ratio compared to their surroundings than sunspots 
in the optical. Therefore, at these wavelengths, when the 
planet transits over an active region, the amount of light 
blocked by the planet becomes greater which produces a 
negative dip in the light curve. There have been efforts in 
the literature to simulate the effects bright spots occul- 
tations and the corresponding dips in light curves have 
on exoplanet transmission spectroscopy (Oshagh et al. 
2014). The light curves in the right panel of Figure 2 
show the effects of the planet transiting over active re¬ 
gions on the transit shape. In the optical, star spots 
are cooler than their surroundings and so appear dark. 
Therefore, when the planet transits over a star spot, the 
fractional loss in light is less and a positive bump is reg¬ 
istered in the light curve (e.g., Rabus et al. 2009). 

The evolution of unocculted active regions on the stel¬ 
lar disk make it difficult to accurately normalize the stel¬ 
lar flux level. Throughout this work we have assumed 
the Sun to be static during a transit observation. Given 
the relatively low activity level, and slow rotation period 
of the Sun, this is a reasonable assumption. It is worthy 
of note however, that for more active stars the corona 
will be evolving rapidly and obtaining a long baseline is 
crucial in accurately normalizing the data and recovering 
an accurate value of Rp/R^. 

3.2. Transit asymmetries 

Figure 3 shows 94 A images of the Sun over two Car¬ 
rington rotations. Under this setup the Sun has rotated 
^ 60° between each image. The top row is for Carring¬ 
ton rotation 2107 (2011 February) when the Sun was in 
activity minimum, and the bottom row is for Carrington 
rotation 2149 (2014 April) when the Sun was more active 
(but still not at activity maximum). The corresponding 
94 A (green) and optical (orange) light curves are shown 
below each image for an impact parameter of 6 = —0.3. 
We find that even for transits separated by a few days 
the 94 A light curves can be highly variable and have 
a different shape when compared to the previous tran¬ 
sit. This is due to the corona being highly structured 
and variable over short timescales. Indeed, for more ac¬ 
tive stars, and for stars with stronger differential rotation 
rates the coronae may be even more dynamic than shown 
in Figure 3 (Pevtsov et al. 2003; Gibb et al. 2014). 

Since high energy transits are sensitive to star light 
from the upper atmosphere of the star including the ex¬ 
tended corona, it may in principle be possible to use tran¬ 
sits at these wavelengths to study the properties of the 
corona. The duration of transit events at these wave¬ 
lengths may be longer than the corresponding optical 
events which are only sensitive to light from the photo¬ 
sphere of the star. In our simulated 94A- 335A light 
curves we see hints of an extended transit when the 
planet has transited over an active region on the stel¬ 
lar limb. We find that the maximum difference in phase 
between the 94 A and optical ingress is ^ 0.005. In prin¬ 
ciple, this effect could be used to study the extent and 
density profile of stellar coronae; however, very high ca¬ 
dence, signal-to-noise transits would be required which 
is currently observationally challenging at these wave¬ 
lengths. 

For WASP-12b, the near-UV (2800 A) ingress and op¬ 
tical ingress timing differed by ^ 0.02 in phase. We found 
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Fig. 3. — Soft X-ray (94 A) images of the Sun over two complete rotations. Over plotted is the orbital path of the simulated planet 
(Rp/R^ = 0.1,6 = —0.3). The corresponding 94 A transit light curves (green) and the simultaneous optical light curve (orange) is shown 
beneath each image. The 94 A light curves reveal how even over short time scales the transits are highly variable due the structure and 
evolution of the corona and rotational modulation. As the planet transits over a bright region of the corona a large dip is registered in the 
94 A light curve but not in the optical light curve. 


no detectable timing difference between our FUV (1600 
& 1700 A) and optical light curves. This therefore sug¬ 
gests that the asymmetry in the light curve of WASP- 
12b is unlikely to be have been caused by stellar activity, 
but indeed could be caused by additional, asymmetri¬ 
cally distributed material surrounding the exosphere of 
the planet. 


3.3. Transit depths 

To further quantify the variability in our simulated 
light curves, we computed the maximum depth AFmax 
for each light curve. Schlawin et al. (2010) provide an es¬ 
timate for the maximum depth of limb-brightened tran¬ 
sits based on a model where the flux from the star is 
primarily concentrated at the limb rather than at disk 
center as is typical for optical wavelengths. Under this 
setup, the maximum expected transit depth (in the ab¬ 
sence of active regions on the stellar surface) occurs at 
second contact of the transit. 

Figure 4 shows a top-down schematic of this setup. 
The length of the arc occulted by the planet is and 
2Rp = R^{1 — cosO) ^ lRi,0‘^ for 6> C 1. The area 



◄-► 

R^-2Rp 


Fig. 4. — Top-down schematic of a limb-brightened transit (not 
drawn to scale) adapted from Schlawin et al. (2010). The image 
shows how the majority of the starlight comes from the limb of the 
star. Under this setup, the deepest part of the transit will be at 
ingress and egress, rather than mid-transit for longer wavelength 
observations. 
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Fig. 5.— Histograms of maximum transit depth, AFmax as a function of wavelength for the 94 A (green), 1700 A (purple), and 4500 A 
(orange) light curves. The left panel is for an impact parameter 6 = 0, the right histogram is for b = —0.3. In the FUV and optical we 
generally recover a consistent value of AFmax (see Table 2). We find a much larger range of maximum transit depths in the soft X-ray 
/ EUV caused by both unocculted active regions creating shallower transits and occulted active regions resulting in deeper transits. The 
results for the other wavelengths are shown in Table 2. 


occulted by the planet can then be approximated by an 
ellipse with semi-major axis and semi-minor axis 

Rp^ implying the occulted area can be expressed as ^4 = 
T:Rp^JRpR^. Finally, the maximum transit depth (in 
normalized units) is then expressed as 

AFbright ^ _A_ ^ 1 

Fbright 27rR2 2[rJ 

For a more detailed geometric description and deriva¬ 
tion of Equation 1 we refer the reader to Schlawin et ah 
(2010). They note that for the same value of Rp/R^, 
transits at wavelengths where the stellar disk is limb- 
brightened will be deeper than those where the star is 
limb-darkened. For a uniformly illuminated stellar disk, 
the maximum transit depth occurs when the entire disk 
of the planet is within the stellar disk (as viewed by the 
observer), and the normalized depth of the transit can 
be expressed as 

AFuniform _ ^^^p _ f ^2) 

^uniform 27rR^ \ R* / 

In the presence of limb-darkening, the mid-transit 
depth will deviate slightly from the depth given in Equa¬ 
tion 2. Equations 1 and 2 provide estimates for Rp/R*; 
however, in practice a full transit model must be fit to 
the light curve to accurately derive Rp/R^ (see Section 
3.5). 

Eigure 5 is a histogram of all values of AEmax com¬ 
puted from all the simulated transit light curves for the 
94 A, 1700 A, and 4500 A light curves. We find that 
AEmax can deviate significantly from the predicted value 
of AEbright given in Equation 1, and that the deepest part 
of the transit may not be at the limb when active regions 
are present on the stellar disk. 

Eor the optical and EUV light curves we find that ac¬ 
tivity has less impact on the depth of the transit and we 
recovered the expected transit depth of 0.01 in ^ 80% 


of the EUV light curves. In the optical, when the planet 
transits over a dark sunspot, the fractional loss in light 
is less and so a relatively localized, positive bump is reg¬ 
istered in the light curve (which is excluded from the fit¬ 
ting). As such, we recovered the expected transit depth 
of 0.01 in all of our optical light curves. Eigure 5 reveals 
how the scatter in AEmax decreases with longer wave¬ 
lengths. 

The largest range in the maximum depth of the transit 
occurs at the shortest wavelengths (94A - 335A). This 
will have significant impact on the ability to reliably re¬ 
cover the radius of the planet (see Section 3.5). The light 
curves where the planet has transited over the equator 
of the Sun (6 = 0) are on average shallower than the 
light curves where the planet transits over the activity 
belts (6 = —0.3). This can be seen in Eigure 6 where 
the average of the ^ 1300 simulated light curves for each 
wavelength in our simulations for both impact parame¬ 
ters is computed. These shallow transits are the direct 
consequence of unocculted bright regions on the stellar 
disk during the transit. In the normalized stellar image 
unocculted bright regions mean the planet occults rel¬ 
atively less star light during the transit, resulting in a 
shallower transit. Eor a more in-depth discussion of un¬ 
occulted dark spots on the depths of transit light curves 
we refer the reader to Section 3 of Pont et al. (2013). 
Note, dark star spots on the transit depth will have the 
opposite effect to the bright regions present at soft X-ray 
and UV wavelengths. 

3.4. Transits affected by activity 

We quantified the percentage of our simulated light 
curves that show signatures of stellar activity (where 
the planet occults an active region) by first simulat¬ 
ing activity-free transits using the code developed by 
Schlawin et al. (2010), itself based on the limb-darkened 
transit code of Mandel & Agol (2002). These transit 
codes assume that the stellar disk is free from active 
regions. We computed a limb-brightened, activity-free. 
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b = 0 b = -0.3 




Fig. 6. — Averaged light curves for our 94 A (green); 171 A (yellow); 304 A (red); 1700 A (purple); 4500 A (orange) simulated light 
curves. The left panel are the averaged light curves for a planet with impact parameter 6 = 0 and the right panel for b = —0.3. The 
limb-brightened light curves in the left panel are shallower than those in the right panel due to the effects of unocculted bright regions on 
the stellar disk. The dashed line is a model limb-brightened, spot free, model transit from Schlawin et al. (2010). 


model transit for a hot Jupiter with Rp/R^ = 0.1 and 
impact parameters 6 = 0 and b = —0.3 (plotted as dashed 
black lines in Figure 6) to match the input to our simula¬ 
tions. For the UV and optical light curves we simulated 
limb-darkened, activity-free transits using the code de¬ 
veloped by Mandel & Agol (2002). 

We then calculated the residuals for our Solar tran¬ 
sits by subtracting the activity-free light curve from each 
simulated light curve to compute the residual flux and 
searched for deviations in the residuals caused by the 
planet occulting an active region. For the X-ray and UV 
transits, when the planet occults a bright active region a 
negative dip is registered in the light curve. We therefore 
searched the residuals for such dips. We conservatively 
classified a dip as being detectable if it’s residual flux 
value was 3x larger than the 1—a uncertainty in the light 
curves of Poppenhaeger et al. (2013), i.e., if the residuals 
contained a value larger than 0.018 x normalized flux. 

When the planet occults an active region in the optical 
a positive bump is registered in the light curve. Guided 
by the measured uncertainties of Sanchis-Ojeda & Winn 
(2011) we classified a bump as being detectable if it’s 
residual flux value was larger than 0.001 x normalized 
flux, making it detectable by Kepler. 

Columns 3 and 5 of Table 2 give the fraction of light 
curves that show detectable signatures of the planet oc¬ 
culting active regions as a function of impact parameter 
(as defined by the criteria above). The fraction of af¬ 
fected light curves is heavily dependent on wavelength 
and geometry of the transit chord with respect to the 
activity belts, with ^ 15% of the 94A - 335A transits 
where 6 = 0 exhibiting detectable signatures of active re¬ 
gion occultation. This fraction increases to ^ 70% when 
the planet transits over the activity belts (6 = —0.3). 
10 — 25% of FUV light curves exhibit detectable dips 
caused by stellar activity depending on the orbital path 
of the planet with respect to the activity belts on the stel¬ 
lar surface. Only ^ 5% of the optical light curves showed 
detectable bumps where the planet had transited over a 
dark sunspot, consistent with the findings of Llama et al. 
(2012). Since these bumps are localized features in the 
light curve they can be accounted for when computing 


Rp/R*- 

More active stars will have a much larger fraction of 
the stellar surface covered with active regions, making 
the likelihood of the planet occulting an active region 
higher. This suggests that our findings are likely to serve 
as lower limits on the fraction of affected light curves for 
active stars. If these dips caused by the planet occulting 
active regions at short wavelengths are not accounted 
for when fitting transit models then an inflated value of 
Rp/R^ will be found. Given the current time sampling 
and sensitivity of high energy transits it is very difficult 
to isolate individual bumps in the light curve. 

3.5. Reeovering Rp/R^ as a funetion of wavelength 

We also recovered the radius ratio, Rp/R^ to inves¬ 
tigate how stellar activity impacts this measurement for 
both the limb-darkened and limb-brightened transits. To 
recover Rp/R^ from the limb-brightened light curves we 
again used the full transit model from Schlawin et al. 
(2010) to produce transit light curves with varying val¬ 
ues of Rp /R^ to find the best fit limb-brightened model 
transit. In the FUV and optical, we used the full transit 
model of Mandel & Agol (2002) with varying values of 
Rp/R* to find the best fit transit. We note that in the 
optical it is relatively straightforward to exclude regions 
of the light curve where the planet has occulted a dark 
star spot; however, at short wavelengths it is non-trivial 
to determine which parts of the light curve have been 
affected by activity and so in the soft X-ray and UV we 
fit the entire light curve. 

Figure 7 shows histograms of the recovered values of 
Rp/R* as a function of wavelength and impact parame¬ 
ter for the 94 A, 1700A, and 4500A transits. Golumns 2 
& 4 of Table 2 give summary statistics of our recovered 
values of Rp/R^ for all the wavelengths. We recovered 
the input value of Rp/R* = 0.1 from all of the optical 
light curves. The mean recovered value of Rp/R^ in the 
FUV was ^ 20% larger than the optical value, caused 
by the active regions appearing as extended, bright re¬ 
gions in the FUV rather than relatively confined, dark 
regions in the optical. Observations and analysis of HD 
209458 have shown that the UV (3000A- 3900A) transit 
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TABLE 2 

Summary of recovered values of Rp/R* and fraction of stellar 

ACTIVITY IMPACTED LIGHT CURVES 



6 = 

0 

6 = 

-0.3 

Wavelength 

Recovered 

Activity 

Recovered 

Activity 

(A) 

Rp/R^ 

Impacted (%) 

Rp/R* 

Impacted (%) 

94 

0.090 ± 0.014 

14 

0.120 ± 0.026 

69 

131 

0.089 ± 0.012 

10 

0.119 ± 0.020 

69 

171 

0.093 ± 0.010 

05 

0.112 ± 0.015 

40 

193 

0.092 ± 0.012 

04 

0.118 ± 0.017 

48 

211 

0.090 ± 0.017 

10 

0.123 ± 0.024 

74 

304 

0.094 ± 0.008 

11 

0.110 ± 0.016 

39 

335 

0.085 ± 0.020 

19 

0.125 ± 0.031 

64 

1600 

0.113 ± 0.007 

10 

0.115 ± 0.009 

24 

1700 

0.118 ± 0.006 

06 

0.117 ± 0.007 

16 

4500 

0.100 ± 0.001 

03 

0.100 ± 0.001 

05 


Note. — Column 1 is the wavelength of the observation. Columns 2 & 3 
give statistics for the simulations with 6 = 0. Columns 4 & 5 show the same 
statistics but for 6 = —0.3. Columns 2 & 4 give the median value of Rp/R* 
and associated variability (derived from the difference between the 90^^ and 
10^^ percentiles of the histograms shown in Figure 7). Columns 3 & 5 give 
the percentage of the light curves where stellar activity caused a detectable 
deviation from the model, activity free light curves (see Section 3.4 for details). 


b = 0.0 b = -0.3 



Fig. 7. — Histograms displaying the recovered values of Rp/R* (see Section 3.5 for more details) for the 94A (green), 1700A (purple), 
and 4500A (orange) simulated light curves. We always recover the input value of Rp/R* =0.1 from the optical light curves. For the 
FUV light curves we find the mean value of Rp/R* to be larger than the optical value due to active regions appearing as extended bright 
regions. For the high energy transits (94A-335A) we find the values of Rp/R* to be more variable. The light curves where the planet did 
not transit over the activity belts (6 = 0) are affected by unocculted active regions which can result in an underestimate in the recovered 
value of Rp/R* of up-to 20%. When the planet transits over the activity belts (6 = —0.3) these light curves are much deeper, resulting in 
an overestimate in the recovered values of Rp/R* of up-to 50%. The blue histogram is the recovered Rp/R* when we rebin our 94A light 
curves to a lower cadence (see Section 4). 


depth is some 0.030% deeper than in the optical (4500A- 
5500A) (Ballester et al. 2007; Sing et al. 2008; Lecave- 
lier Des Etangs et al. 2008). Lecavelier Des Etangs et al. 
(2008) attribute this additional absorption to Rayleigh 
scattering in the atmosphere of the planet. Eor such 
small changes in transit depth it is therefore imperative 
that the effects of unocculted active regions (which can 
also change the transit depth) are accounted for in UV 
transit light curve analysis. 

There was a much wider spread in the recovered ra¬ 
dius values at the shortest wavelengths. The radius of 
the planet was underestimated by up-to 25% in the light 
curves where the planet did not transit over the activity 
belts (5 = 0). This underestimate is due to a shallower 
transit caused by the presence of unocculted bright re¬ 


gions on the stellar disk. When the planet did transit 
over the activity belts (5 = —0.3), the recovered value 
of Rp/R^ was found to be over-estimated by up-to 50%, 
caused by the planet transiting over bright regions re¬ 
sulting in a deeper transit. This can be seen in Eigure 
8 where we have plotted the best fit X-ray light curve 
(dashed blue line) for two light curves. Again, for more 
active stars these values will serve as lower limits. 

4. IMPLICATIONS FOR X-RAY TRANSITS 

X-ray (lA - 63A) transit observations of HD 189733b 
have revealed a transit depth up-to 4.8 times deeper than 
the optical transit depth which has been attributed to the 
planet having an extended atmosphere (Poppenhaeger 
et al. 2013). These X-ray light curves have a mid-transit 































Stellar activity impact on X-ray and UV transits 


9 


X 

D 


T3 

a> 

"fU 

E 


2011 January 17 


1.00 

0.98 

0.96 

0.94 



0.98 1.00 1.02 

Phase 


2013 December 26 



0.98 1.00 1.02 

Phase 


Fig. 8. — The optical light curve (orange) and 94 A (green) for nearly-quiescent (left) and activity impacted (right) cases. The best 
fit transit light curve to the X-ray light curve is shown as the dashed blue line. In the lower-left of each panel is the corresponding 94 A 
image. The black points are the result of rebinning the 94 A light curve and adding a similar level of noise to the levels in the observations 
of Poppenhaeger et al. (2013). We find that the maximum depth of the 94 A low-cadence light curves is typically shallower than for the 
high-cadence light curves; however, fewer data points make it more challenging to distinguish between an increased transit being caused 
by stellar activity, or by an extended planetary atmosphere. 


depth shallower than the depth in the optical. In our sim¬ 
ulated short wavelength (94 A - 335 A) transits, where 
Rp/R* is kept constant, only those light curves unaf¬ 
fected by stellar activity have mid-transit depths shal¬ 
lower than the optical depth (see the right panel in Fig¬ 
ure 2). This suggests that either HD 189733b may indeed 
have an extended atmosphere in X-rays, or that during 
those observations the planet was transiting over X-ray 
bright regions that covered a large fraction of the transit 
path. Both of these scenarios would result in an X-ray 
mid-transit depth deeper than the corresponding optical 
transit depth. This second scenario is not unreasonable 
given that HD 189733 is a more active star than the Sun. 
Indeed, Zeeman-Doppler Imaging of HD 189733 has re¬ 
vealed the star has a large-scale magnetic field strength 
ten times stronger than that of the Sun (Fares et al. 
2010 ). 

More X-ray transits will therefore help determine 
which scenario is correct. We find that for large active 
regions, a bump in the optical light curve coincides with 
large dips in the high energy light curves, suggesting that 
optical light curves could potentially be used to rule out 
stellar activity effects on deep X-ray light curves. There¬ 
fore, simultaneous optical light curves will also aide the 
interpretation. 

Due to the difficulty in obtaining X-ray observations, 
the light curves of Poppenhaeger et al. (2013) are at a 
lower cadence than those simulated here. In our high- 
cadence observations the signatures of the planet tran¬ 
siting over an active region can be resolved as localized 
dips in the light curve; however, this may not always 
be the case for lower cadence observations. The black 
points in Figure 8 show the effects of rebinning two 94A 
light curves to the same cadence of the X-ray transits of 
Poppenhaeger et al. (2013) and adding similar sized er¬ 
ror bars. These rebinned light curves make it much more 
difficult to detect the signatures of the planet transiting 


over active regions. The right panel of Figure 8 shows 
how the localized dips in the high cadence light curve 
caused by the planet occulting an active region are less 
defined at lower cadence. If the planet occults larger (or 
more) active regions than in this case then the depth of 
the X-ray transit will be consistently lower than the op¬ 
tical transit depth. This may lead to the false conclusion 
that the planet has an inflated atmosphere in X-rays. 

To quantify the effects of using lower cadence observa¬ 
tions we repeated the analysis carried out in Section 3.5 
and derived the value of Rp/R^ using the low cadence 94 
A light curves. The dark blue histogram in Figure 5 gives 
the results for this analysis. We find a reasonable agree¬ 
ment between the recovered values of Rp/R* for the light 
curves when the planet is not occulting the activity belts 
(6 = 0). When the planet is occulting the activity belts 
(6 = —0.3), the low cadence light curves produce smaller 
values of Rp/R^ when compared to the full cadence ob¬ 
servations. This is caused by the active regions being 
localized regions, and so integrating for longer during 
the observations reduces their impact on the light curve; 
however, the disadvantage of having fewer data points is 
that it is less obvious whether the increased depth in the 
light curve is caused by the planet or stellar activity. 

5. CONCLUSIONS 

In this paper we used resolved images of the Solar disk 
to simulate transit light curves at multiple wavelengths 
ranging from soft X-rays/EUV to optical. By simulating 
a transiting planet with two impact parameters we inves¬ 
tigated the importance of both occulted and unocculted 
active regions on transit light curves. 

At short wavelengths, up-to 70% of the simulated tran¬ 
sit light curves where the planet occulted the activity 
belts on the Sun showed signatures of stellar activity 
that would result in the derived value of Rp/R^ deviating 
from the real value. This value decreased to ~ 20% for 
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the light curves when the planet’s orbital path did not oc¬ 
cult the activity belts. At longer wavelengths, the effects 
of stellar activity were not as prevalent, with 20% of our 
simulated FUV light curves where the planet occulted 
the activity belts showing signatures of stellar activity. 
In the optical, we found that for stars with activity lev¬ 
els similar to the Sun, only 5% of light curves will have 
bumps caused by the planet transiting over star spots 
agreeing with the findings of Llama et al. (2012). 

We also investigated the effects stellar activity would 
have on recovering the radius ratio of the planet and the 
host star, Rp/R^. If the planet transits over the activity 
belts on the star, then the mean value of the planet’s 
radius as derived from soft X-ray and EUV light curves 
can be overestimated by 20%; however, during times of 
high activity the radius can be overestimated by up-to 
50%. When the planet did not occult the activity belts, 
the effects of unocculted active regions resulted in the 
mean value of the planet’s radius at short wavelengths 
being underestimated by 10%. Again, during periods of 
high activity this underestimate increased to ^ 25%. 

The effects were similar in the FUV for which the ra¬ 
dius of the planet could be over-estimated by ^ 15%. In 
the optical we recovered the input value of Rp/R^ = 0.1 
in all our simulated light curves since bumps in the light 
curve can be accounted for when deriving Rp/R^. 

We compared our light curves to other transit obser¬ 
vations at short wavelengths. We found that transits 
at short wavelengths showed hints of an extended dura¬ 
tion than the corresponding optical transits. This is due 
to these observations being sensitive to light from both 
the photosphere and extended corona, suggesting high- 


energy transits could be used to study the density and 
structure of stellar coronae. Unlike the limb-brightened 
transits, the FUV (1600 A & 1700 A) light curves exhib¬ 
ited no signatures of an early ingress when compared to 
the optical. This suggests that the observed early-ingress 
in the near-UV (2800 A) transit of WASP-12 is unlikely 
to be caused by stellar activity alone; rather, it strength¬ 
ens the idea that WASP-12b has additional absorbing 
material surrounding the exosphere of the planet. 

By rebinning our high-cadence soft X-ray (94 A) light 
curves to a similar cadence as the X-ray data of HD 
189733b taken by Poppenhaeger et al. (2013), we demon¬ 
strated that it is difficult to distinguish between a deeper 
X-ray transit being caused by stellar activity, or by an 
extended planetary atmosphere. Higher time sampling 
will make identifying parts of the light curve where the 
planet has occulted an active region easier. Simultane¬ 
ous optical data may also aide in revealing whether the 
planet is transiting over active regions and help distin¬ 
guish between these two scenarios. 

We conclude that it is very important that when at¬ 
tempting to determine planetary properties using tran¬ 
sits at high energy wavelengths, multiple observations 
are obtained with long baselines, and that simultaneous 
optical data (both during the transit and out-of-transit) 
may help decouple planetary properties from stellar vari- 
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